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Abstract In this lecture we aim to present a formalism based on Faddeev-like 
equations for describing nuclear three-cluster reactions that include elastic, trans- 
fer and breakup channels. Two different techniques based on momentum-space and 
configuration-space representations are explained in detail. An important new fea- 
ture of these methods is the possibility to account for the repulsive Coulomb inter- 
action between two of the three clusters in all channels. Comparison with previous 
calculations based on approximate methods used in nuclear reaction theory is also 
discussed. 



1 Introduction 

Nuclear collision experiments, performed at ion accelerators, are a very powerful 
tool to study nuclear properties at low and intermediate energies. In order to inter- 
pret accumulated experimental data appropriate theoretical methods are necessary 
enabling the simultaneous description of the available elastic, rearrangement and 
breakup reactions. 

Regardless of its importance, the theoretical description of quantum-mechanical 
collisions turns out to be one of the most complex and slowly advancing problems in 
theoretical physics. If during the last decade accurate solutions for the nuclear bound 
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State problem became available, full solution of the scattering problem (containing 
elastic, rearrangement and breakup channels) remains limited to the three-body case. 

The main difficulty is related to the fact that, unlike the bound state wave func- 
tions, scattering wave functions are not localized. In configuration space one is 
obliged to solve multidimensional differential equations with extremely complex 
boundary conditions; by formulating the quantum-mechanical scattering problem 
in momentum space one has to deal with non-trivial singularities in the kernel of 
multivariable integral equations. 

A rigorous mathematical formulation of the quantum mechanical three-body 
problem in the framework of non relativistic dynamics has been introduced by Fad- 
deev in the early sixties [19], in the context of the three-nucleon system with short 
range interactions. In momentum space these equations might be slightly modified 
by formulating them in terms of three-particle transition operators that are smoother 
functions compared to the system wave functions. Such a modification was proposed 
by Alt, Grassberger, and Sandhas [1] (AGS). 

Solutions of the AGS equations with short range interactions were readily ob- 
tained in the early seventies. As large computers became available progress fol- 
lowed leading, by the end eighties, to fully converged solutions of these equations 
for neutron-deuteron (n-d) elastic scattering and breakup using realistic short range 
nucleon-nucleon (N-N) interactions. Nevertheless the inclusion of the long range 
Coulomb force in momentum space calculations of proton-deuteron (p-d) elastic 
scattering and breakup with the same numerical reliability as calculations with short 
range interactions alone, only become possible in the last decade. 

Significant progress has been achieved [16, 17] by developing the screening and 
renormalization procedure for the Coulomb interaction in momentum space using a 
smooth but at the same time sufficiently rapid screening. This technique permitted to 
extend the calculations to the systems of three-particles with arbitrary masses above 
the breakup threshold [11, 18]. 

However it has taken some time to formulate the appropriate boundary conditions 
in configuration space for the three-body problem [32, 33, 36] and even longer to 
reformulate the original Faddeev equations to allow the incorporation of long-range 
Coulomb like interactions [34, 35]. Rigorous solution of the three-body problem 
with short range interactions has been achieved just after these theoretical develop- 
ments, both below and above breakup threshold. On the other hand the numerical 
solution for the three-body problem including charged particles above the three- 
particle breakup threshold has been achieved only recently. First it has been done 
by using approximate Merkuriev boundary conditions in configuration space [27]. 
Nevertheless this approach proved to be a rather complex task numerically, remain- 
ing unexplored beyond the p-d scattering case, but not yet for the p-d breakup. 

Finally, very recently configuration space method based on complex scaling have 
been developed and applied for p-d scattering [31]. This method allows to treat the 
scattering problem using very simple boundary conditions, equivalent to the ones 
employed to solve the bound-state problem. 
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The aim of this lecture is to present these two recently developed techniques, 
namely the momentum-space method based on screening and renormalization as 
well as the configuration-space complex scaling method. This lecture is structured 
as follows: the first part serves to introduce theoretical formalisms for momentum 
space and configuration space calculations; in the second part we present some se- 
lected calculations with an aim to test the performance and validity of the two pre- 
sented methods. 



2 Momentum-space description of three-particle scattering 

We describe the scattering process in a system of three-particles interacting via pair- 
wise short-range potentials Va, 01 = 1 , 2, 3; we use the odd-man-out notation, that is, 
vi is the potential between particles 2 and 3. In the framework of nonrelativistic 
quantum mechanics the center-of-mass (cm.) and the internal motion can be sepa- 
rated by introducing Jacobi momenta 

_ rriykp-mpky 

Pa — ; , (.i) 

mp + niy 

nia (kft + ky) - (mfi + my)ka 
Qa = : , (2) 

with (ajSy) being cyclic permutations of (123); ka and nia are the individual particle 
momenta and masses, respectively. The cm. motion is free and in the following we 
consider only the internal motion; the corresponding kinetic energy operator is Hq 
while the full Hamiltonian is 

3 

H^Ho+Y.Va. (3) 
a=l 



2.1 Alt, Grassberger, and Sandhas equations 

We consider the particle a scattering from the pair a that is bound with energy £«. 
The initial channel state \ba({a) is the product of the bound state wave function \ba) 
for the pair a and a plane wave with the relative particle-pair a momentum ({a, the 
dependence on the discrete quantum numbers is suppressed in our notation. \ba(ia) 
is the eigenstate of the corresponding channel Hamiltonian Ha = Hq + Vu with the 
energy eigenvalue E = Ea + q^j^/lMa where Ma is the particle-pair a reduced mass. 
The final channel state is the particle-pair state in the same or different configuration 
I^/bQ/j) in the case of elastic and rearrangement scattering or, in the case of breakup, 
it is the state of three free particles |pyqy) with the same energy E = /9y/2/iy + 
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c^yjlMy and pair yreduced mass \iy, any set of Jacobi momenta can be used equally 
well for the breakup state. 

The stationary scattering states [42, 22] corresponding to the above channel states 
are eigenstates of the full Hamiltonian; they are obtained from the channel states 
using the full resolvent G = (£ + /O — //)^\ i.e., 

|fcaqa>'+* = /0G|/7aq„), (4) 
|paqa)'+* = /OGlpaqa). (5) 

The full resolvent G may be decomposed into the channel resolvents G^~(E^- 
iO — Hp)^^ and/or free resolvent Go = (E + iO — Hq)^^ as 

G = Gp+GpvpG, (6) 

with p =0,1,2,3 and v^g = Ly=i ^pyVy where 5py^ 1 — Spy. Furthermore, the chan- 
nel resolvents 

Gp = Go + GoTpGo, (7) 
can be related to the corresponding two-particle transition operators 

Tp=vp+vpGoTp, (8) 

embedded into three-particle Hilbert space. Using these definitions Eqs. (4) and (5) 
can be written as triads of Lippmann-Schwinger equations 

|/'aqa)<+'= 5pa\ba(la)+GpVp\bccqa)^+\ (9) 
|Paqa)<+'= il+GoTp)\paqa)+GpVp\paqa)^+\ (10) 

with a being fixed and j3 = 1,2,3; they are necessary and sufficient to define the 
states l^'aqa)'-^' and Ipaqa)'^' uniquely. However, in scattering problems it may 
be more convenient to work with the multichannel transition operators Upa defined 
such that their on-shell elements yield scattering amplitudes, i.e., 

Upa\baqa) ^vplbaqa)'-^^- (H) 

Our calculations are based on the AGS version [1] of three-particle scattering theory. 
In accordance with Eq. (11) it defines the multichannel transition operators Upa by 
the decomposition of the full resolvent G into channel and/or free resolvents as 

SpaGa + GpUpaGa- (12) 

The multichannel transition operators Up^ with fixed a and j8 = 1 , 2, 3 are solutions 
of three coupled integral equations 

3 

Upa = ^PaGQ +Y^5pyTyGQiUya. (13) 

y=l 
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The transition matrix Uoa to final states with three free particles can be obtained 
from the solutions of Eq. (13) by quadrature, i.e., 

3 

Uoa^G^' + Y^TyGaUya. (14) 

y=i 

The on-shell matrix elements {bpt(p\Upa\baqa) are ampUtudes (up to a factor) 
for elastic (/3 = a) and rearrangement ()3 7^ a) scattering. For example, the differ- 
ential cross section for the a + (/37) j3 + (ya) reaction in the cm. system is given 
by 

do q' 

= (2;r)4M„M^-^|(Z7^q'^|f/p„|/7„q„)|2. (15) 

The cross section for the breakup is determined by the on-shell matrix elements 
(Pylyl^Oal^ceQa)- Thus, in the AGS framework all elastic, rearrangement, and 
breakup reactions are calculated on the same footing. 

Finally we note that the AGS equations can be extended to include also the three- 
body forces as done in Ref. [12]. 



2.2 Inclusion of the Coulomb interaction 

The Coulomb potential wc, due to its long range, does not satisfy the mathemat- 
ical properties required for the formulation of standard scattering theory as given 
in the previous subsection for short-range interactions v^. However, in nature the 
Coulomb potential is always screened at large distances. The comparison of the 
data from typical nuclear physics experiments and theoretical predictions with full 
Coulomb is meaningful only if the full and screened Coulomb become physically 
indistinguishable. This was proved in Refs. [44, 43] where the screening and renor- 
malization method for the scattering of two charged particles was proposed. We base 
our treatment of the Coulomb interaction on that idea. 

Although we use momentum-space framework, we first choose the screened 
Coulomb potential in configuration-space representation as 

w^(r)=wc(r)e-('-/'^)", (16) 

and then transform it to momentum-space. Here R is the screening radius and « con- 
trols the smoothness of the screening. The standard scattering theory is formally 
applicable to the screened Coulomb potential wr, i.e., the Lippmann-Schwinger 
equation yields the two-particle transition matrix 

tR^WR + WRgQtR, (17) 

where go is the two-particle free resolvent. It was proven in Ref. [44] that in the 
limit of infinite screening radius R the on-shell screened Coulomb transition matrix 
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(screened Coulomb scattering amplitude) (p'|fi;|p) with 75' = p, renormalized by an 
infinitely oscillating phase factor z]f^{p) — e^"l''^^P\ approaches the full Coulomb 
amplitude (p'|fc|p) in general as a distribution. The convergence in the sense of 
distributions is sufficient for the description of physical observables in a real exper- 
iment where the incoming beam is not a plane wave but wave packet and therefore 
the cross section is determined not directly by the scattering amplitude but by the 
outgoing wave packet, i.e., by the scattering amplitude averaged over the initial state 
physical wave packet. In practical calculations [2, 16] this averaging is carried out 
implicitly, replacing the renormalized screened Coulomb amplitude in the R ^ °o 
limit by the full one, i.e., 

limz„'(/,)(p'|f«|p)^(p'|fc|p). (18) 

K — 

Since Zg ' (p) is only a phase factor, the above relations indeed demonstrate that the 
physical observables become insensitive to screening provided it takes place at suf- 
ficiently large distances R and, in the R ^ °° limit, coincide with the corresponding 

_ 1 

quantities referring to the full Coulomb. Furthermore, renormalization by - (/?, ) 
in the R^ °o limit relates also the screened and full Coulomb wave functions [23], 
i.e., 

lim (1 +8otR)\p)zr-ip) = \wt\p)). (19) 

K — >oo 

The screening and renormalization method based on the above relations can be 
extended to more complicated systems, albeit with some limitations. We consider 
the system of three-particles with charges Za of equal sign interacting via pairwise 
strong short-range and screened Coulomb potentials Va + WaR with a being 1, 2, 
or 3. The corresponding two-particle transition matrices are calculated with the full 
channel interaction 

7^*' = {va + WaR) + {va + WaR)Gi)T^'' , (20) 

and the multichannel transition operators U^^ for elastic and rearrangement scat- 
tering are solutions of the AGS equation 

3 

^/3a ^ ^PuGq^ + ^ ^PyTy GoUya (21) 

r=i 

all operators depend parametrically on the Coulomb screening radius R. 

In order to isolate the screened Coulomb contributions to the transition amplitude 
that diverge in the infinite R limit we introduce an auxiliary screened Coulomb po- 
tential W^'^^' between the particle a and the center of mass (cm.) of the remaining 
pair. The same screening function has to be used for both Coulomb potentials Wqr 
and W^^'. The corresponding transition matrix 

r„T = + wxclf'r^r, (22) 
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with Ga = {E + iO — Hq — Va ^ Wur)^^ is a two-body-like operator and therefore 
its on-shell and half-shell behavior in the limit R^°°is given by Eqs. (18) and (19). 
As derived in Ref. [16], the three -particle transition operators may be decomposed 
as 

<j = Sp^T^J^- + [1 + r^TGjf [1 + ^« '^«"] (23) 

= 5^„r^r + (<j-5/3ar^r)- (24) 

where the auxihaiy operator U^^ is of short range when calculated between on- 

shell screened Coulomb states. Thus, the three -particle transition operator U^J has 

a long-range part SpaT^^ whereas the remainder Up^ — ^fiaJaR' is a short-range 
operator that is externally distorted due to the screened Coulomb waves generated 

(R\ 

by [1 + Ga TaR ']- On-shell, both parts do not have a proper limit as oo but 
the limit exists after renormalization by an appropriate phase factor, yielding the 
transition amplitude for full Coulomb 

{bpc(p\uf^\ba^a) = haMp\^ac\ba^a) 

+ hm^[Z^7 {4p ) {bpvip I (Ijfl - ^p^TaR ) l^aqa)Z„i • (25) 

The first term on the right-hand side of Eq. (25) is known analytically [44]; it cor- 
responds to the particle-pair a full Coulomb transition amplitude that results from 
the impUcit renormalization of T^^- according to Eq. (18). The °° limit for the 

remaining part {ujf^ — SpaT^^ ') of the multichannel transition matrix is performed 
numerically; due to the short-range nature of this term the convergence with the 
increasing screening radius R is fast and the limit is reached with sufficient accu- 
racy at finite R; furthermore, it can be calculated using the partial-wave expansion. 
We emphasize that Eq. (25) is by no means an approximation since it is based on 
the obviously exact identity (24) where the /? — > oo limit for each term exists and is 
calculated separately. 

The renormalization factor for R ^ °ois a diverging phase factor 

Za«(^a)-e-2'*««('?°), (26) 

where 'PaRi^a), though independent of the particle-pair relative angular momentum 
la in the infinite R limit, may be realized by 

^afl (?a ) = < (^a ) - {qa), (27) 

with the diverging screened Coulomb phase shift '7/"^(^a) corresponding to stan- 
dard boundary conditions and the proper Coulomb one o"^{qa) referring to the 
logarithmically distorted proper Coulomb boundary conditions. For the screened 
Coulomb potential of Eq. (16) the infinite R limit of 'PaRiqa) is known analytically. 
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<i'aR{qa) = ■Jfra{qa)[\n{2qaR) ^ C /n], (28) 

where Csa 0.5772 156649 is the Euler number and J^^(^a) = Oie.m.Za Y^y^yaZyMal qa 
is the Coulomb parameter with ag.m. ~ 1/137. The form of the renormalization 
phase 'PaRiqa) to be used in the actual calculations with finite screening radii R 
is not unique, but the converged results show independence of the chosen form of 
'^aRiqa)- 

For breakup reactions we follow a similar strategy. However, the proper three- 
body Coulomb wave function and its relation to the three-body screened Coulomb 
wave function is, in general, unknown. This prevents the application of the screening 
and renormalization method to the reactions involving three free charged particles 
(nucleons or nuclei) in the final state. However, in the system of two charged parti- 
cles and a neutral one with Zp = 0, the final-state Coulomb distortion becomes again 
a two-body problem with the screened Coulomb transition matrix 

TpR = WpR + WpRGoTpR. (29) 

This makes the channel p, corresponding to the correlated pair of charged particles, 
the most convenient choice for the description of the final breakup state. As shown 
in Ref. [17], the AGS breakup operator 

t^da = Go' + t7^f>Gof/<^\ (30) 

y=l 

can be decomposed as 

Ul^a = (1 +7p«Go)C'(l (31) 

where the reduced operator L^q^^ (Z) calculated between screened Coulomb distorted 
initial and final states is of finite range. In the full breakup operator UqJ (Z) the ex- 

(R) 

ternal distortions show up in screened Coulomb waves generated by (1 + Ga T^r') 
in the initial state and by (1 + TprGq) in the final state; both wave functions do not 
have proper limits as /? ^ oo. Therefore the full breakup transition amplitude in the 
case of the unscreened Coulomb potential is obtained via the renormalization of the 
on-shell breakup transition matrix U^J in the infinite R limit 

(P^Qp It^Oa l^aq«> = jim [Zpi {p'p){v'pi{p |^aqa>zj (^a)], (32) 

where p'p is the relative momentum between the charged particles in the final state, 
Qp the corresponding particle-pair relative momentum, and 

ZpRip'p) = e-2''^p(''p)li"(2/^«)-c/«l, (33) 
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the final-state renormalization factor with the Coulomb parameter Kp{pp) for the 
pair p. The limit in Eq. (32) has to be performed numerically, but, due to the short- 
range nature of the breakup operator, the convergence with increasing screening 
radius R is fast and the limit is reached with sufficient accuracy at finite R. Thus, to 
include the Coulomb interaction via the screening and renormalization method one 
only needs to solve standard scattering theory equations. 



2.3 Practical realization 

We calculate the short-range part of the elastic, rearrangement, and breakup scat- 
tering amplitudes (25) and (32) by solving standard scattering equations (21), (22), 
and (30) with a finite Coulomb screening radius R. We work in the momentum-space 
partial-wave basis [10], i.e., we use three sets 

\paqaVa) = \paqa{la{[La{si3Sy)Sa]IaSa}Ka)JM) with {a,p,Y) being cyclic per- 
mutations of (1,2,3). Here Sa is the spin of particle a. La and la are the orbital an- 
gular momenta associated with pa and qa respectively, whereas Sa, la, and Ka are 
intermediate angular momenta that are coupled to a total angular momentum J with 
projection M. All discrete quantum numbers are abbreviated by Va - The integration 
over the momentum variables is discretized using Gaussian quadrature rules thereby 
converting a system of integral equations for each J and parity FI = (— )^a+/a jjjj-Q 
a very large system of linear algebraic equations. Due to the huge dimension those 
linear systems cannot be solved directly. Instead we expand the AGS transition op- 
erators (21) into the corresponding Neumann series 

<j = ^aG,' + f BpyTl^^Sya + i: a^^rf >Go SraT^'^^B^a + • • • , (34) 
7=1 7=1 (7=1 

that are summed up by the iterative Pade method [5]; it yields an accurate solution 
of Eq. (21) even when the Neumann series (34) diverges. Each two-particle tran- 

( R) 

sition operator Ty is evaluated in its proper basis IpygyVy), thus, transformations 
between all three bases are needed. The calculation of the involved overlap functions 
IPalaVa) follows closely the calculation of three-nucleon permutation op- 
erators discussed in Refs. [10, 22]. A special treatment [5, 10] is needed for the in- 
tegrable singularities arising from the pair bound state poles in Ty^' and from Gq. 
Furthermore, we have to make sure that R is large enough to achieve (after renormal- 
ization) the i?-independence of the results up to a desired accuracy. However, those 
R values are larger than the range of the nuclear interaction resulting in a slower 
convergence of the partial-wave expansion. As we found in Ref. [16], the practi- 
cal success of the screening and renormalization method depends very much on the 
choice of the screening function, in our case on the power n in Eq. (16). We want to 
ensure that the screened Coulomb potential wr approximates well the true Coulomb 
one wc for distances r < R and simultaneously vanishes rapidly for r > R, provid- 
ing a comparatively fast convergence of the partial-wave expansion. As shown in 



10 



A. Deltuva, A. C. Fonseca, and R. Lazauskas 



Fig. 1 Differential cross 
section and deuteron vec- 
tor analyzing power (Ti i of 
the ad elastic scattering at 
4.81 MeV deuteron lab en- 
ergy as functions of the cm. 
scattering angle. Convergence 
with the screening radius R 
used to calculate the short- 
range part of the amplitudes 
is studied: R = 5 fm (dotted 
curves), R = 10 fm (dash- 
dotted curves), and R = 15 fm 
(solid curves). Results without 
Coulomb are given by dashed 
curves. The experimental data 
are from Refs. [4, 30]. 




60 120 180 

©cm. (deg) 



Ref. [16], this is not the case for simple exponential screening (n = 1) whereas the 
sharp cutoff (n ^ oo) yields slow oscillating convergence with the screening radius 
R. However, we found that values of 3 < n < 8 provide a sufficiently smooth and 
rapid screening around r ^ R. The screening functions for different n values are 
compared in Ref. [16] together with the results demonstrating the superiority of our 
optimal choice: using 3 < « < 8 the convergence with the screening radius R, at 
which the short range part of the amplitudes was calculated, is fast enough such 
that the convergence of the partial-wave expansion, though being slower than for 
the nuclear interaction alone, can be achieved and there is no need to work in a 
plane-wave basis. Here we use n = 4 and show in Figs, land 2 few examples for the 
/?-convergence of the a-deuteron scattering observables calculated in a three-body 
model («,/?,«); the nuclear interaction is taken from Ref. [11]. The convergence 
with R is impressively fast for both a-deuteron elastic scattering and breakup. In 
addition we note that the Coulomb effect is very large and clearly improves the 
description of the experimental data, especially for the differential cross section in 
a-deuteron breakup reaction. This is due to the shift of the ap P-wave resonance 
position when the a p Coulomb repulsion is included that leads to the corresponding 
changes in the structure of the observables. 

In addition to the internal reliability criterion of the screening and renormal- 
ization method — the convergence with R — we note that our results for proton- 
deuteron elastic scattering [15] agree well over a broad energy range with those 
of Ref. [28] obtained from the variational configuration-space solution of the three- 
nucleon Schrodinger equation with unscreened Coulomb potential and imposing the 
proper Coulomb boundary conditions explicitly. 
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Fig. 2 Fivefold differen- 
tial cross section of the 
ad breakup reaction at 
15 MeV a lab energy for 
several combinations of a 
and proton scattering an- 
gles as function of the final- 
state energy variable S with 
dS = {dEl+dEf,y/^. Con- 
vergence with the screening 
radius R is studied: R= 10 fm 
(dotted curves), S = 15 fm 
(dash-dotted curves), and 
R = 20 fm (solid curves). 
Results without Coulomb are 
given by dashed curves. The 
experimental data are from 
Ref. [29]. 
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3 Configuration space 

In contrast to the momentum-space representation, the Coulomb interaction has a 
trivial expression in configuration space and thus may seem to be easier to han- 
dle. However the major obstacle for configuration-space treatment of the scattering 
problem is related with the complexity of the wave function asymptotic structure, 
which strongly complicates once three-particle breakup is available. Although for 
short range interactions the analytical behavior of the breakup asymptote of the con- 
figuration space wave function is well established, this is not a case once long range 
interactions (like Coulomb) are present. Therefore a method which enables the scat- 
tering problem to be solved without explicit use of the wave function asymptotic 
form is of great importance. The complex scaling method has been proposed [39, 6] 
and successfully applied to calculate the resonance positions [37] by using bound 
state boundary conditions. As has been demonstrated recently this method can be ex- 
tended also for the scattering problem [7, 46]. We demonstrate here that this method 
may be also successfully applied to solve three-particle scattering problems which 
include the long-range Coulomb interaction together with short range optical poten- 
tials. 



3.1 Faddeev- Merkuriev equations 

Like in the momentum space formalism described above Jacobi coordinates are 
also used in configuration space to separate the center of mass of the three-particle 
system. One has three equivalent sets of three-particle Jacobi coordinates 
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' 2mpmy 
(mp + my)m 



[ry-rp), (35) 



/ 2mp{mp+my) mprp+myTy 

ya = \ 7 ; ; (^a ; ), 

y (w^a + >np + niyjin mp + my 

here Va and ma are individual particle position vectors and masses, respectively. The 
choice of a mass scale m is arbitrary. The three-particle problem is formulated here 
using Faddeev-Merkuriev (FM) equations [34]: 



(£ - //o - ^ vv') V/a = (va + Wa){Wa + + Wy). 
3 

[E -Ho-Y, »'i)Wi3 = (vp + Wp){Wa + WI3 + Wr). (36) 

3 

{E-Hq- W'i)\j/y = {Vy + Wy){\jfa + Wl3 + ¥r)^ 

K=l 

where the Coulomb interaction is split in two parts (short and long range), Wa = 
+ i^a' ^y means of some arbitrary cut-off function Xa{xa,ya)' 

w'aixa.ya) = Waixa)Xaixa.ya) w'^ixa^ya) = Wa{xa)[l - Xa{xa.ya)] (37) 

This cut-off function intends to shift the full Coulomb interaction in the term 
if Xa is small, whereas the term acquires the full Coulomb interaction if Xa 
becomes large and ya < Xa- The practical choice of function Xa{xa,ya) has been 
proposed in [34]: 

, (38) 

with free parameters xo,yo having size comparable with the charge radii of the re- 
spective binary systems; the value of parameter jU must be larger than 1 and is usu- 
ally set jj. Ril.ln such a way the so-called Faddeev amplitude xj/a intends to acquire 
full asymptotic behavior of the binary a — (jSy) channels, i.e: 

ra(Xa,ya = 5^,aVaH^a)^a''" (ya) + YfjaiA^cc.yaWtM^t-""'iya) 

Ja 

+ A„(x„,y„)0-"(p), (39) 

where the hyperradius is p = \/^a +3'a- An expression (Pa {xa)<j)'a''" (ya) repre- 
sents the incoming wave for particle a on pair (fiy) in the bound state /«, with 
(p'a{xa) representing the normalized wave function of bound state /«. This wave 
function is a solution of the (£ — //q — — \'a — Wa'"' ) two-body Hamiltonian. The 
^a'''""{ya) and <P""'{pa) represent outgoing waves for binary and three-particle 
breakup channels respectively. In the asymptote, one has the following behavior: 
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(Paixa ^ °°) °<: exp(-fc,„JCa), 
^a'""'{ya oo)oc exp{iqi„ya), (40) 
^Cip ^^)'^exp{iKp), (41) 

with kj^ — y^—e.^m representing momentum of 2-body bound state with a nega- 
tive binding energy e.^ ; qt^ = ^{E — e.^ )m is relative scattering momentum for the 

a — {^y) binary channel, whereas K = \/mE is a three-particle breakup momentum 
(three-particle breakup is possible only if energy value E is positive). 

When considering particle's a scattering on the bound state ia of the pair (jSy), 
it is convenient to separate readily incoming wave ^f'a''" = W'a i'^a)^a '"' {y a)^ by 
introducing: 

Wi?'""' = W'a"-Wa" i^aWa'" (ya), (42) 
Then Faddeev-Merkuriev equations might be rewritten in a so-called driven form: 





-Ho- 


k:=1 


= (v'ct - 










{E 


-Ho~ 


k:=1 


= i^r 










{E 


-Ho~ 


k:=1 













In this expression index of the incoming state has been omitted in all Faddeev 
component expressions ij/'^ and i//^"'. 



3.2 Complex scaling 

Next step is to perform the complex scaling operations i.e. scale all the distances x 
and y by a constant complex factor e'^ , so that both Re{e'^) and Im{e'^) are posi- 
tive (angle 6 must be chosen in the first quartet in order to satisfy this condition). 
The complex scaling operation, in particular, implies that the analytical continua- 
tion of the interaction potentials is performed: x'aixae'^) and Waixae'^). Therefore 
the complex scaling method may be used only if these potentials are analytic. It 
is easy to see that the solutions of the complex scaled equations coincide with the 
ones obtained without complex scaling but to which the complex scaling operation 
is applied: [\i/{xa,ya)f^ = r(x„e''',3;„e''^). 

Namely, it is easy to demonstrate that all the outgoing wave functions of eq.(41) 
becomes exponentially bound after the complex scaling operation: 
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[(Pa {xa ^ °°)]'^'^ °= exp(-fc/„XaCOS0), 

r / out i*-^ 

<l>a' ' {ya^°°) exp(-^,„yasin0), (44) 
[0,7 (p ^-)]^^ - expi-Kpsine). 
Nevertheless an incoming wave diverges in ya after the complex scaling: 

<Pa'""'{ya -^°°) °^ exp{+qi„ya sin 0). (45) 

However these terms appear only on the right hand sides of the driven Faddeev- 
Merkuriev equation (44) being pre-multiplied with the potential terms and under 
certain conditions they may vanish outside of some finite (resolution) domain Xa G 
OjX™"] and ya G [Oj^™^"]. Let us consider the long range behavior of the term 

{vp+ Wp ) xjf'a ■ Since the interaction terms vp and are of short range, the only 

region the former term might not converge is along yp axis in {xp,yp ) plane, i.e. for 

xp <^yp - On the other hand Xai^pjyp) ~ yJmy/{my + mp)y^M/ {my + ma)yp and 

ya{iip,yp) ~ ^ mp/{my + mp ) yjma/{my + ma)yp under condition ^yp. Then 
one has: 



oc exp -ki^ 



nijM 



[niy + mp)(my + ma) 



yp cos 9 + fy,„ 



{my + mp){mj + ma)' 
(46) 

This term becomes bound to finite domain in {xp,yp) plane, if condition: 



-vb sin 9 



tan0 < 



j MyM 

mainp 





B 


a 




E + 


B 1 

• a 1 



(47) 



is satisfied. This implies that for rather large scattering energies £, above the break- 
up threshold, one is obliged to use rather small complex scaling parameter Q values. 

^f'a'"', in principle, is not exponentially bound 



The term [^^^ 



Wa- 



after the complex scaling. It represents the higher order corrections to the residual 
Coulomb interaction between particle a and bound pair ()3 7). These corrections are 
weak 0(1/3'^) and might be neglected by suppressing this term close to the border 
of the resolution domain. Alternative possibility might be to use incoming wave 
functions, which account not only for the bare a — ()3 7) Coulomb interaction but 
also takes into account higher order polarization corrections. 



Extraction of the scattering observables is realized by employing Greens theo- 
rem. One might demonstrate that strong interaction amplitude for a — (/3 7) collision 
is: 
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Hja 



cs 



(wry iya+wa-w^"'r'm 



with i^'i^f = 



cs 



(48) 



being the total wave function 



of the three-body system. In the last expression the term containing product of two 
incoming waves is slowest to converge. Even stronger constraint than eq.(47) should 
be implied on complex scaling angle in order to make this term integrable on the 
finite domain. Nevertheless this term contains only the product of two-body wave 
functions and might be evaluated without using complex scaling prior to three-body 
solution. Then the appropriate form of the integral (48) to be used becomes: 



cs 



Ja-in 
a 



CS 



m 



{¥'an*i^a + wa-wr)¥t'"d'xid'y 



e^'^d^Xid'^yi 



(49) 



4 Application to three-body nuclear reactions 

The two methods presented in sections 2 and 3 were first applied to the proton- 
deuteron elastic scattering and breakup [16, 17, 12, 31]. The three-nucleon system 
is the only nuclear three-particle system that may be considered realistic in the sense 
that the interactions are given by high precision potentials valid over a broad energy 
range. Nevertheless, in the same way one considers the nucleon as a single particle 
by neglecting its inner quark structure, in a further approximation one can consider 
a cluster of nucleons (composite nucleus) to be a single particle that interacts with 
other nucleons or nuclei via effective potentials whose parameters are determined 
from the two-body data. A classical example is the a particle, a tightly bound four- 
nucleon cluster. As shown in Figs. 1 and 2 and in Ref. [11], the description of the 
{a,p,n) three-particle system with real potentials is quite successful at low energies 
but becomes less reliable with increasing energy where the inner structure of the a 
particle cannot be neglected anymore. At higher energies the nucleon-nucleus or 
nucleus-nucleus interactions are modeled by optical potentials (OP) that provide 
quite an accurate description of the considered two-body system in a given narrow 
energy range; these potentials are complex to account for the inelastic excitations 
not explicitly included in the model space. The methods based on Faddeev/AGS 
equations can be applied also in this case, however, the potentials within the pairs 
that are bound in the initial or final channel must remain real. The comparison of the 
two methods based on the AGS and FM equations will be performed in section 4. 1 
for such an interaction model with OP. 

In the past the description of three-body-like nuclear reactions involved a number 
of approximate methods that have been developed. Well-known examples are the 
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distorted-wave Born approximation (DWBA), various adiabatic approaches [26], 
and continuum-discretized coupled-channels (CDCC) method [3]. Compared to 
them the present methods based on exact Faddeev or AGS equations, being more 
technically and numerically involved, have some disadvantages. Namely, their ap- 
phcation in the present technical realization is so far limited to a system made of 
two nucleons and one heavier cluster The reason is that the interaction between two 
heavier cluster involves very many angular momentum states and the partial-wave 
convergence cannot be achieved. The comparison between traditional nuclear reac- 
tion approaches and momentum-space Faddeev/AGS methods for various neutron + 
proton + nucleus systems are summarized in section 4.2. 

On the other hand, the Faddeev and AGS methods may be more flexible with 
respect to dynamic input and thereby allows to test novel aspects of the nuclear 
interaction not accessible with the traditional approaches. Few examples will be 
presented in section 4.3. 



4.1 Numerical comparison of AGS and FM methods 

As an example we consider the n + p+^-^C system. For the n-p interaction we use 
a realistic AVI 8 model [47] that accurately reproduces the available two-nucleon 
scattering data and deuteron binding energy. To study not only the d + '^C but also 
p+ '^C scattering and transfer reactions we use a n-'^C potential that is real in 
the ^Pi partial wave and supports the ground state of '^C with 4.946 MeV binding 
energy; the parameters are taken from Ref. [38]. In all other partial waves we use 
the n-^^C optical potential from Ref. [45] taken at half the deuteron energy in the 
d + '^C channel. The p-^^C optical potential is also taken from Ref. [45], however, 
at the proton energy in the p + ^^C channel. We admit that, depending on the reaction 
of interest, other choices of energies for OP may be more appropriate, however, the 
aim of the present study is comparison of the methods and not the description of the 
experimental data although the latter are also included in the plots. 

We consider d + '^C scattering at 30 MeV deuteron lab energy and p + ^^C scat- 
tering at 30.6 MeV proton lab energy; they correspond to the same energy in cm. 
system. First we perform calculations by neglecting the p-^^C Coulomb repulsion. 
One observes a perfect agreement between the AGS and FM methods. Indeed, the 
calculated S-matrix elements in each three-particle channel considered (calculations 
have been performed for total three -particle angular momentum states up to 7 = 13) 
agree within three digits. Scattering observables converge quite slowly with J as 
different angular momentum state contributions cancel each other at large angles. 
Nevertheless, the results of the two methods are practically indistinguishable as 
demonstrated in Fig. 3 for d + '^C elastic scattering and transfer to p + '^C. 

Next we perform the full calculation including the p-'^C Coulomb repulsion; 
we note that inside the nucleus the Coulomb potential is taken as the one of a uni- 
formly charged sphere [11]. Once again we obtain good agreement between the 
AGS and FM methods. However, this time small variations up to the order of 1% 
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are observed when analyzing separate 5-matrix elements, mostly in high angular 
momentum states. This leads to small differences in some scattering observables, 
e.g., differential cross sections for d + '^C elastic scattering (at large angles where 
the differential cross section is very small) and for the deuteron stripping reaction 
d+'^^C p+'^^C shown in Fig. 4. The p + elastic scattering observables pre- 
sented in Fig. 5 converge faster with J. As a consequence, the results of the two 
calculations are indistinguishable for the /? + '^C elastic cross section and only tiny 
differences can be seen for the proton analyzing power at large angles. In any case, 
the agreement between the AGS and FM methods exceeds both the accuracy of the 
data and the existing discrepancies between theoretical predictions and experimental 
data. 



Fig. 3 Comparison of 
momentum- (solid curves) 
and configuration-space 
(dashed-dotted curves) re- 
sults for the deuteron-'^C 
scattering at 30 MeV deuteron 
lab energy. Differential cross 
sections for elastic scatter- 
ing and stripping are shown 
neglecting the Coulomb inter- 
action. 
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4.2 Comparison with traditional nuclear reaction approaches 

The method based on the momentum-space AGS equations has already been used 
to test the accuracy of the traditional nuclear reaction approaches; limitations of 
their validity in energy and kinematic range have been estalished. The distorted- 
wave impulse approximation for breakup of a one-neutron halo nucleus ^'Be on 
a proton target has been tested in Ref. [9] while the adiabatic-wave approxima- 
tion for the deuteron stripping and pickup reactions "Be(;?,c/)"'Be, '^C(c/,/7)'^C, 
and 48Ca(^f,p)49c a in Ref. [38]. However, one of the most sophisticated traditional 
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Fig. 4 Comparison of 
momentum- (solid curves) 
and configuration-space 
(dashed-dotted curves) re- 
sults for the deuteron-'^C 
scattering at 30 MeV deuteron 
lab energy. Differential cross 
sections for elastic scattering 
and stripping are shown, the 
former in ratio to the Ruther- 
ford cross section daRjdQ. 
The experimental data are 
fromRefs. [41,40]. 
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Fig. 5 Comparison of 
momentum- (solid curves) 
and configuration-space 
(dashed-dotted curves) results 
for the proton- '^^C elastic 
scattering at 30.6 MeV pro- 
ton lab energy. Differential 
cross section divided by the 
Rutherford cross section and 
proton analyzing power are 
shown. The experimental data 
are from Ref. [25]. 
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approaches is the CDCC method [3]. A detailed comparison between CDCC and 
AGS results is performed in Ref. [18]. The agreement is good for deuteron-'^C 
and deuteron-^^Ni elastic scattering and breakup. In these cases nucleon-nucleus 
interactions were given by optical potentials; thus, there was no transfer reaction. 
A different situation takes place in proton-"Be scattering where '^Be nucleus is 
assumed to be the bound state of a '''Be core plus a neutron. In this case, where 
the transfer channel d + '"Be is open, the CDCC approach lacks accuracy as shown 
in Ref. [18]. The semi-inclusive differential cross section for the breakup reaction 
p + "Be p + n + '"Be was calculated also using two CDCC versions where the 
full scattering wave function was expanded into the eigenstates of either the n + '"Be 
(CDCC-BU) or the p + n (CDCC-TR) pair. Neither of them agrees well with AGS 
over the whole angular regime as shown in Fig. 6. It turns out that, depending on 
the '"Be scattering angle, the semi-inclusive breakup cross section is dominated 
by different mechanisms: at small angles it is the proton-neutron quasifree scatter- 
ing whereas at intermediate and large angles it is the neutron-'"Be D-wave reso- 
nance. However, a proper treatment of proton-neutron interaction in CDCC-BU and 
of neutron- '"Be interaction in CDCC-TR is very hard to achieve since the wave 
function expansion uses eigenstates of a different pair No such problem exists in 
the AGS method that uses simultaneously three sets of basis states and each pair is 
treated in its proper basis. 




Fig. 6 Semi-inclusive differential cross section for the breakup reaction p+ "Be — )■ '"Be 
at lab energy of 38.4 MeV/nucleon. Results obtained with AGS and CDCC methods are compared. 
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4.3 Beyond standard dynamic models 

The standard nucleon-nucleus optical potentials employed in three-body calcula- 
tions have central and, eventually, spin-orbit parts that are local. This local approx- 
imation yields a tremendous simplification in the practical realization of DWBA, 
CDCC and other traditional approaches that are based on configuration-space rep- 
resentations where the use of nonlocal optical potentials was never attempted. How- 
ever, nonlocal optical potentials do not yield any serious technical difficulties in the 
momentum-space representation. Thus, they can be included quite easily in the AGS 
framework employed by us. 

There are very few nonlocal parametrizations of the optical potentials available. 
We take the one from Refs. [20, 21] defined in the configuration space as 

vy{r\r)^H,{x)[V,{y) + iW,.{y)]+2SyLyH,{x)V,{y), (50) 

with X ~ |r' — r| and y = |r' + r|/2. The central part has real volume and imag- 
inary surface parts, whereas the spin-orbit part is real; all of them are expressed 
in the standard way by Woods-Saxon functions. Some of their strength parameters 
were readjusted in Ref. [13] to improve the description of the experimental nucleon- 
nucleus scattering data. The range of the nonlocality is determined by the functions 
Hi{x) = (7rj3?)^^/^exp(— x^/j3?) with the parameters j3,- being of the order of 1 fm. 

A detailed study of nonlocal optical potentials in three-body reactions involving 
stable as well as weakly bound nuclei, ranging from '"Be to ^"Ca, is carried out in 
Ref. [13]. In order to isolate the nonlocality effect we also performed calculations 
with a local optical potential that provides approximately equivalent description of 
the nucleon-nucleus scattering at the considered energy. The nonlocality effect turns 
out to be very small in the elastic proton scattering from the bound neutron-nucleus 
system and of moderate size in the deuteron-nucleus scattering. However, the ef- 
fect of nonlocal proton-nucleus optical potential becomes significant in deuteron 
stripping and pickup reactions {d,p) and {p,d); in most cases it considerably im- 
proves agreement with the experimental data. Examples for {d,p) reactions leading 
to ground and excited states of the stable nucleus '^O and one-neutron halo nu- 
cleus '^C are presented in Figs. 7 and 8. We note that in these transfer reactions the 
proton-nucleus potential is taken at proton lab energy in the proton channel while 
the neutron-nucleus potential has to be real in order to support the respective bound 
states. 

Another extension beyond the standard dynamic models includes the AGS method 
using energy-dependent optical potentials Although such calculations don't corre- 
spond to a rigorous Hamiltonian theory, they may shed some light on the short- 
comings of the traditional nuclear interaction models. A detailed discussion of the 
calculations with energy-dependent optical potentials is given in Ref. [14]. 
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Fig. 7 Differential cross 
section for {d, p) reaction on 
^''O at 36 MeV deuteron lab 
energy leading to "O nucleus 
in the ground state 5/2+ 
(top) and first excited state 
1/2+ (bottom). Predictions 
of nonlocal (solid curve) and 
local (dashed curve) optical 
potentials (OP) are compared 
with the experimental data 
from Ref. [8]. 
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Fig. 8 Differential cross 
section for (d,p) reaction 
on '*C at 14 MeV deuteron 
lab energy leading to one- 
neutron halo nucleus '^C in 
the ground state 1/2+ (top) 
and first excited state 5/2+ 
(bottom). Curves as in Fig. 7 
and the experimental data are 
from Ref [24]. 
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5 Summary 

We have presented the results of three-body Faddeev-type calculations for systems 
of three particles, two of which are charged, interacting through short-range nuclear 
plus the long-range Coulomb potentials. Realistic applications of three-body the- 
ory to three-cluster nuclear reactions — such as scattering of deuterons on a nuclear 
target or one-neutron halo nucleus impinging on a proton target — only became pos- 
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sible to address in recent years when a reliable and practical momentum-space treat- 
ment of the Coulomb interaction has been developed. After the extensive and very 
complete study of p-d elastic scattering and breakup, the natural extension of these 
calculations was the application to complex reactions such as li-^He, p-^^O, ' 'Be-/5, 
J-^^Ni and many others using a realistic interaction such as AVI 8 between nucle- 
ons, and optical potentials chosen at the appropriate energy for the nucleon-nucleus 
interactions. The advantage of three-body calculations vis-a-vis traditional approx- 
imate reaction methods is that elastic, transfer, and breakup channels are treated on 
the same footing once the interaction Hamiltonian has been chosen. Another advan- 
tage of the three-body Faddeev-AGS approach is the possibility to include nonlocal 
optical potentials instead of local ones as commonly used in the standard nuclear 
reaction methods; as demonstrated, this leads to an improvement in the description 
of transfer reactions in a very consistent way across different energies and mass 
numbers for the core nucleus. 

Although most three-body calculations have been performed in momentum space 
over a broad range of nuclei from ^He to ^**Ni and have encompassed studies of cross 
sections and polarizations for elastic, transfer, charge exchange, and breakup reac- 
tions, coordinate space calculations above breakup threshold are coming to age us- 
ing the complex scaling method. We have demonstrated here that both calculations 
agree to within a few percent for all the reactions we have calculated. This is a very 
promising development that may bring new light to the study of nuclear reactions 
given that the reduction of the many-body problem to an effective three-body one 
may be better implemented and understood by the community in coordinate space 
rather than in momentum space. On the other hand, compared to DWBA, adiabatic 
approaches, or CDCC, the Faddeev-type three-body methods are computationally 
more demanding and require greater technical expertise rendering them less attrac- 
tive to analyze the data. Nevertheless, when benchmark calculations have been per- 
formed comparing the Faddeev-AGS results with those obtained using CDCC or 
adiabatic approaches, some discrepancies were found in transfer and breakup cross 
sections depending on the specific kinematic conditions. Therefore the Faddeev- 
AGS approach is imminent in order to calibrate and validate approximate nuclear 
reaction methods wherever a comparison is possible. 
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